Robust Conditional Independence maps of single-voxel Magnetic Resonance Spectra to elucidate associations between brain tumours and metabolites.

The aim of the paper is two-fold. First, we show that structure finding with the PC algorithm can be inherently unstable and requires further operational constraints in order to consistently obtain models that are faithful to the data. We propose a methodology to stabilise the structure finding process, minimising both false positive and false negative error rates. This is demonstrated with synthetic data. Second, to apply the proposed structure finding methodology to a data set comprising single-voxel Magnetic Resonance Spectra of normal brain and three classes of brain tumours, to elucidate the associations between brain tumour types and a range of observed metabolites that are known to be relevant for their characterisation. The data set is bootstrapped in order to maximise the robustness of feature selection for nominated target variables. Specifically, Conditional Independence maps (CI-maps) built from the data and their derived Bayesian networks have been used. A Directed Acyclic Graph (DAG) is built from CI-maps, being a major challenge the minimization of errors in the graph structure. This work presents empirical evidence on how to reduce false positive errors via the False Discovery Rate, and how to identify appropriate parameter settings to improve the False Negative Reduction. In addition, several node ordering policies are investigated that transform the graph into a DAG. The obtained results show that ordering nodes by strength of mutual information can recover a representative DAG in a reasonable time, although a more accurate graph can be recovered using a random order of samples at the expense of increasing the computation time

Tracking Therapy Response in Glioblastoma Using 1D Convolutional Neural Networks

Background: Glioblastoma (GB) is a malignant brain tumour that is challenging to treat, often relapsing even after aggressive therapy. Evaluating therapy response relies on magnetic resonance imaging (MRI) following the Response Assessment in Neuro-Oncology (RANO) criteria. However, early assessment is hindered by phenomena such as pseudoprogression and pseudoresponse. Magnetic resonance spectroscopy (MRS/MRSI) provides metabolomics information but is underutilised due to a lack of familiarity and standardisation. Methods: This study explores the potential of spectroscopic imaging (MRSI) in combination with several machine learning approaches, including one-dimensional convolutional neural networks (1D-CNNs), to improve therapy response assessment. Preclinical GB (GL261-bearing mice) were studied for method optimisation and validation. Results: The proposed 1D-CNN models successfully identify different regions of tumours sampled by MRSI, i.e., normal brain (N), control/unresponsive tumour (T), and tumour responding to treatment (R). Class activation maps using Grad-CAM enabled the study of the key areas relevant to the models, providing model explainability. The generated colour-coded maps showing the N, T and R regions were highly accurate (according to Dice scores) when compared against ground truth and outperformed our previous method. Conclusions: The proposed methodology may provide new and better opportunities for therapy response assessment, potentially providing earlier hints of tumour relapsing stages

Anti-tumour immune response in GL261 glioblastoma generated by Temozolamide Immune-Enhancing Metronomic Schedule monitored with MRSI-based nosological images

Glioblastomas (GB) are brain tumours with poor prognosis even after aggressive therapy. Improvements in both therapeutic and follow-up strategies are urgently needed. In previous work we described an oscillatory pattern of response to Temozolomide (TMZ) using a standard administration protocol, detected through MRSI-based machine learning approaches. In the present work, we have introduced the Immune-Enhancing Metronomic Schedule (IMS) with an every 6-days TMZ administration at 60mg/kg, and investigated the consistence of such oscillatory behaviour. A total of n=17 GL261 GB tumour-bearing C57BL/6j mice were studied with MRI/MRSI every 2 days, and the oscillatory behaviour (6.2±1.5 days period from the TMZ administration day) was confirmed during response. Furthermore, IMS-TMZ produced significant improvement in mice survival (22.5±3.0 days for controls vs 135.8±78.2 for TMZ-treated), outperforming standard TMZ treatment. Histopathological correlation was investigated in selected tumour samples (n=6) analyzing control and responding fields. Significant differences were found for CD3+ cells (lymphocytes, 3.3±2.5 vs 4.8±2.9 respectively) and Iba-1 immunostained area (microglia/macrophages, 16.8±9.7% and 21.9±11.4% respectively). Unexpectedly, during IMS-TMZ treatment, tumours from some mice (n=6) fully regressed and remained undetectable without further treatment for one month. These animals were considered “cured” and a GL261 re-challenge experiment performed, with no tumour reappearance in 5 out of 6 cases. Heterogeneous therapy response outcomes were detected in tumour-bearing mice, and a selected group was investigated (n=3 non-responders, n=6 relapsing tumours, n=3 controls). PD-L1 content was found ca. 3-fold increased in the relapsing group when comparing with control and non-responding groups, suggesting that increased lymphocyte inhibition could be associated to IMS-TMZ failure. Overall, data suggest that host immune response has a relevant paper in therapy response/escape in GL261 tumours under IMS-TMZ therapy. This is associated to changes in the metabolomics pattern, oscillating every 6 days, in agreement with immune cycle length, which is being sampled by MRSI-derived nosological images

A Novel Semi-Supervised Methodology for Extracting Tumor Type-Specific MRS Sources in Human Brain Data

BackgroundThe clinical investigation of human brain tumors often starts with a non-invasive imaging study, providing information about the tumor extent and location, but little insight into the biochemistry of the analyzed tissue. Magnetic Resonance Spectroscopy can complement imaging by supplying a metabolic fingerprint of the tissue. This study analyzes single-voxel magnetic resonance spectra, which represent signal information in the frequency domain. Given that a single voxel may contain a heterogeneous mix of tissues, signal source identification is a relevant challenge for the problem of tumor type classification from the spectroscopic signal.Methodology/Principal FindingsNon-negative matrix factorization techniques have recently shown their potential for the identification of meaningful sources from brain tissue spectroscopy data. In this study, we use a convex variant of these methods that is capable of handling negatively-valued data and generating sources that can be interpreted as tumor class prototypes. A novel approach to convex non-negative matrix factorization is proposed, in which prior knowledge about class information is utilized in model optimization. Class-specific information is integrated into this semi-supervised process by setting the metric of a latent variable space where the matrix factorization is carried out. The reported experimental study comprises 196 cases from different tumor types drawn from two international, multi-center databases. The results indicate that the proposed approach outperforms a purely unsupervised process by achieving near perfect correlation of the extracted sources with the mean spectra of the tumor types. It also improves tissue type classification.Conclusions/SignificanceWe show that source extraction by unsupervised matrix factorization benefits from the integration of the available class information, so operating in a semi-supervised learning manner, for discriminative source identification and brain tumor labeling from single-voxel spectroscopy data. We are confident that the proposed methodology has wider applicability for biomedical signal processing

A multi-centre, web-accessible and quality control-checked database of in vivo MR spectra of brain tumour patients

No description supplie